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SUMMARY 

The laminar boundary -layer flow about a circular cone at large 
angles of attack to a supersonic stream has been analyzed in the plane 
of symmetry by a method applicable in general to the flow about conical 
bodies . 

At the bottom of the cone, velocity profiles were obtained showing 
the expected tendency of the boundary layer to become thinner on the 
under side of the cone as the angle of attack is increased. 

At the top of the cone, the analysis failed to yield unique solu- 
tions, except for small angles of attack. Beyond a certain critical 
angle of attack, boundary-layer flow does not exist in the plane of sym- 
metry, thus indicating separation. This critical angle is presented as 
a f unction of Mach number and cone vertex angle. 


INTRODUCTION 

The supersonic aerodynamics of pointed bodies has considerable cur- 
rent interest in connection with the design of aircraft and missile fuse- 
lages . An important feature of the flow about such bodies is the behavior 
of the boundary layer and, in particular, the flow separation which may 
occur along the t low-pressure side of the body due to angle of attack. 

The present report will consider the development of the laminar boundary 
layer on the surface of a right circular cone at an angle of attack to a 
supersonic stream (see fig. l) . The conical configuration may be con- 
sidered an idealization of the nose portion of a supersonic aircraft 
fuselage . 

Outside a thin boundary layer on a cone, the nonviscous supersonic 
flow (upon which the boundary layer itself depends) is "conical" in the 
sense that physical quantities (such as velocity and pressure) are con- 
stant along any ray proceeding from the cone apex. The description of 
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this outer flow, contained in references 1 to 3 and elucidated in refer- 
ence 4, is considered adequate for the purposes of this report, hut sub- 
ject to restrictions which will be discussed subsequently. 

In figure 2 is shown qualitatively the circumferential pressure dis- 
tribution on the cone surface predicted for various angles of attack 
(see reference 4). These pressure distributions depend only on the char- 
acter of the nonviscous flow beyond the boundary layer, on the assumption 
that the boundary layer is extremely thin* When the angle of attack is 
very small, the pressure decreases monotonically from the bottom of the 
cone around to the top. For larger angles of attack there appears a 
region near the top of the cone wherein the pressure gradient reverses 
and the pressure increases toward the top. As the angle of attack is 
further increased, this region becomes greater in extent. 

As a consequence of the conical nature of the nonviscous flow, it 
is shown in references 5 and 6 that the laminar boundary layer has para- 
bolic similarity along generators of the cone; that is, velocity, pres- 
sure,, and density inside the boundary layer are constant along any para- 
bola (i see fig. l) drawn in any one meridional plane (plane passing through 
the body axis). Of course, circumferential variation of these quantities 
is to be expected when the cone is at angle of attack. 

In reference 7 the effect of angle of attack on the laminar boundary 
layer is analyzed, in the limit of very small angle of attack, with the 
result that the boundary layer tends to be thicker on the top of the cone 
than on the bottom (fig. 3(a)). This is to be expected since the fluid 
near the base of the boundary layer has low inertia, is therefore 
inclined to follow the direction of the circumferential pressure gra- 
dient more closely than is the outer flow, and thus tends to drain away 
from beneath the cone and accumulate near the top. No separation is 
encountered because, for small angle of attack, the pressure gradient is 
always favorable (fig. 2(a)). 

For larger angles of attack, when the pressure gradient reverses 
direction near the top of the cone (fig. 2(b)), experiment indicates the 
formation of boundary-layer "lobes" (fig. 3(b)). When the angle of 
attack is further increased, the lobe pattern finally breaks away from 
the body to form a vortex street (fig. 3(c)). 

Recently (reference 8), at the Lewis laboratory a brief experiment 
was carried out in which a total -head probe was placed near the surface 
at the top of a cone and pointed toward the cone apex (fig. 4(a)). The 
cone was mounted in a supersonic wind tunnel, and the probe was used to 
measure the total head in the boundary layer at a fixed height above 
the surface as the angle of attack was varied by rotating the cone in 
the meridional plane containing the probe. Figure 4(b) shows the 
result of this test. The decrease in indicated total head as the angle 
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of attack was increased from negative to positive values may be inter- 
preted to mean an increase in boundary-layer thickness at the top of the 
cone. Beyond a certain angle of attack, this tendency reverses, and the 
boundary layer apparently becomes thinner as the angle of attack is 
increased. This is a clear indication of the tendency to form lobes, as 
illustrated in figure 3(b) . 

In the present report of research conducted at the NACA Lewis lab- 
oratory, the laminar boundary layer in the meridional plane of symmetry 
of the flow is analyzed for large angles of attack in order to provide 
velocity profiles on the bottom of the cone and to provide a certain 
degree of insight into the question of separation on the top. 


SOLUTION OF BOUNDARY-LAYER EQUATIONS IN PLANE OF SYMMETRY 
Boundary-Layer Equations in Plane of Symmetry 


In reference 7, it is shown that the dimensionless laminar boundary - 
layer equations for supersonic flow over a circular cone are 


1_ p'((p) _2_ 

f + 30 p S + 30 Stp 


2 2 

f XX - 30 SX f \cp + 3(sx) 2 + 2f XXX = 0 ( la ) 


f + 2_Pikl, + 2 

30 p g 30 g <P| 


2 2^2 p'(<p) 

g XX " 30 g X g X(p- 3 g X f X 30 p + 2 g XXX~ 0 

(lb) 

I 2 + -I 2 (10 


T + (f\) 2 + (g\) 2 = T-l + 


(Id) 

Equations (la) and (lb) are momentum equations, equation (lc) is an 
energy balance, and equation (id) is the equation of state. A com- 
plete list of symbols is provided in appendix A. The functions 
f(X,<p) and g(X,<p) are related to the two-component vector poten- 
tial discussed in reference 6 and axe defined according to the 
relat ions 

u = f x 
w = gx 

in a manner such as to satisfy the continuity, equation identically. 

The coordinate X has been formed as follows: 



r-i 

2r 


pT 
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• v, r- 

The coordinate <p is the- angle between the vertical plane of symmetry 
of the flow and any meridional plane of the body (fig. l) . Equations (l) 
and (3) imply that parabolic similarity of the Blasius type exits in 
meridional planes. As pointed out in references 5 and 6, this conclusion 
of parabolic similarity applies for the boundary layer on any smooth 
conically symmetric body in supersonic flow (for example, a cone of ellip- 
tic cross section) . 

In reference 7, all quantities are made dimensionless by referring 
them to the properties of the nonviscous flow at the outer edge of the 
boundary layer when the cone is at zero angle of attack. In the present 
report it will be convenient to use a different reference condition (sub- 
script r) which will be defined subsequently. The following quantities 
on the left are to be identified with the dimensionless groups on the 
right : 

u,w 

P 

T 

P 

x,y ~ p r u r x/Cp r , p r u r y/Cp r J 

where the constant C arises from the assumption of the temperature - 
viscosity relation of Chapman and Rubes in (reference 9): 


u/ur* v / u r ^ 

p/p r 
2c p T / u r 2 
p/p ^ 2 


(*) 





(5a) 


with C being defined as follows, in order to match e'quat ion (5a) to the 
Sutherland formula at the cone surface (denoted by subscript w) : 


C = 



1 + S/T r 
^w/^r + 5/T r 


The quantity S may be taken as equal to (216° R)2c /u 


(5b) 
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The following additional physical assumptions are embodied in equa- 
tions (l): 

(a) A thin boundary layer across which the static pressure is 
constant 

(b) Prandtl number of 1 and constant ratio of specific heats y 

(c) No heat transfer through the surface 

From equation (lc), since the case of Prandtl number of 1 and no 
heat transfer is considered, T w in equation (5b) may be taken equal to 
the dimensionless stream stagnation temperature . 

The boundary conditions on the functions f(X,<p) and g(X, (p) are: 

At the outer edge of the boundary layer, the u and w velocity compo- 
nents should take on the corresponding nonviscous values 

= ui(<p) (6a) 

■ w x (<p) (6b) 

At the cone surface, the u and w velocities should vanish 

fx(°^) = sx(0,<p) = o (6c) 

and the normal, velocity v should vanish. It is shown in reference 6 

that this last requirement is met if 

f(0,<p) = g(0,(p) = 0 (6d) 

Equations (l), involving two independent variables, would be quite 
difficult to solve in general. However, a certain amount of information 
can be obtained by restricting consideration to the plane of symmetry, 
thus yielding a tractable set of ordinary equations involving X as the 
only independent variable. 

In the plane of symmetry cp = 0 ,jc, w = g^= 0. Because u = is 
even about the plane of symmetry and may be expected to be regular there, 
f\q> = 0. The pressure and the density are also even, and therefore 

P'(<p) vanishes at <p = 0,rt. Thus, in the plane of symmetry, equation (la) 
reduces to the following equation: 



(7a) 



6 


MCA TN 2844 


Every term in equation (lb) vanishes at the plane of symmetry; and, there- 
fore, in order to obtain a meaningful equation, it is necessary first to 
differentiate equation (lb) with respect to (p and then drop terms which 
vanish at <P = 0,it . This procedure yields the following result : 


f + 30 Scp/SXXcp " 30 (§X<p) " 3 SXcp^X “ 30 


2 P ”(<P) 


2 §XXX<p 


Equation (lc) becomes 


T + (f x )^ = T-l + u-l 2 (7c) 

Equations (7) may be considered a set of ordinary differential equa- 
tions for the functions f(X,0) and g (p (X,0), or f(X,jt) and g<p(X,rt), 

depending on whether the solution is required at the bottom or the top of 
the cone. According to equations (2), the result f(X,0 or it) may be 
differentiated with respect to X to give the profile of meridional vel- 
ocity u in the plane of symmetry. The form assumed by the circumfer- 
ential velocity profile w as q> •+ 0 or jt is given by or rt) 

in the sense that, at a small angular distance dc p away from the plane of 
symmetry , w * d< p. 

The boundary conditions (equations (6)) become, in the plane of 
symmetry, 

f x (~,0 or 3t) = Up(0 or jt) 

Scp^O or Jt) = w-L (0 or it) S(8) 

f x(0,0 or rt) = g(px(0,0 or si) - f(0,0 or it) = £^>(0,0 or it) = Oj 

In view of the first of equations (8) and of equation (7c), it is con- 
venient to specify the reference condition (subscript r) to be that 
existing at the outer edge of the boundary layer, for the particular 
angle of attack under consideration, evaluated at either <p = 0 or 
(p = rt, depending on whether the analysis pertains to the bottom or top 
of the cone. 

Because the pressure is assumed constant across the boundary layer, 
equations (id) and (7c) and the assumption of constant pressure across 
the boundary layer (p = Pp) provide that 


= i + i 
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From conventions (4), Tj_ = 
definitions axe made: 



For convenience, the following 


Scp(X,0 

or it) = It k>|r(X) 

(lOa) 

k = 

2 

3 e ^C 0 or *0 

(10b) 


Equation (9) and definitions (10) are introduced into equations (7a) and 
(7b), and a value of p’’( 0 or it) is assigned consistent with the non- 
viscous equations at the outer edge of the boundary layer (p''(0 or it) 
may conveniently be obtained from equations (lb) and (6) by setting 

ScpXX = s<t>\\x = o when x=~). 


The following pair of simultaneous ordinary differential equations 
then result: 


(f + kf )f ' * + 2f " ' =0 

f 

(lla) 

(f + ta|r)lr" + 2i|r f " - k(i') 2 - | *'f' + (k + |j\l + i-jl 

\ v, X 

1 

i^h n \ 

ii 

o 


; (nb) 

and boundary conditions (8) become 


f '(-) = ♦'(•) = 1 

(12a) 

f'(0) = +'(0) = 0 

(l2b) 

f(o) = t(0) = 0 

(12c) 


Two parameters appear: k, which depends essentially on angle of 

attack, and Tj_, which is essentially dependent on Mach number. If the 
angle of attack (and hence k) is zero or nearly zero, equations (ll) 
become precisely those considered in reference 7 and may be solved quite 
readily, since equation (lib) becomes linear and the solution of equa- 
tion (lla) is well known as the Blasius function. When k differs sub- 
stantially from zero (moderate or large angle of attack), equations (ll) 
are both nonlinear and the solutions are interdependent. For any par- 
ticular case, when only the stream Mach number, cone vertex angle, and 
angle of attack are specified, the parameters k and Tj must be 
obtained by recourse to a theory of the outer nonviscous flow. 


Outer Nonviscous Flow 

In references 1 to 3, the results of a theory of nonviscous super- 
sonic flow about circular cones at angle of attack are tabulated. The 
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case of zero angle of attack (reference l) is solved exactly in the sense 
that no assumption of small vertex angle is made. The equations are then 
expanded in powers of angle of attack with the use of the zero -angle -of - 
attack solution as the first approximation. Terms linear in angle of 
attack are presented in reference 2, and terms proportional to the square 
of angle of attack are presented in reference 3. Reference 4 clarifies 
the application of the theory to the computation of flow conditions at 
the cone surface . 


There are two objections to the use of this theory in the present 
application: 


(1) Neglect of terms in the expansion beyond that involving the 
square of angle of attack may lead to an insufficiently accurate repre- 
sentation of the flow at the large angles of attack which are of inter- 
est. Unfortunately ; no nonviscous theory is available that treats the 
effect of angle of attack with greater precision. In reference 4, a com- 
parison of the theoretical and experimental pressure distributions is pre- 
sented for a cone of semi vertex angle of 10°, at a Mach number of 2, and 
of an angle of attack of 12.2°. The agreement shown is very good, 
especially since the angle of attack is sufficiently large that the pres- 
sure distribution is of the type shown in figure 2(b). 

(2) In reference 10, Ferri points out that the method of expansion 
used in references 2 and 3 is improper near the cone surface and leads 

to an erroneous form of the entropy distribution around the cone. There- 
fore, the theory cannot be applied if the vertex angle, the angle of 
attack, or the Mach number is so large that the flow may not be considered 
essentially isentropic. In reference 7, an argument is presented to the 
effect that in the limit of infinitesimal angle of attack the presence of 
a boundary layer ensures that the error in entropy distribution is of no 
consequence even for large cone vertex angles. That argument in no way 
applies to the present analysis because the angles of attack considered 
are not infinitesimal. For the .purposes of this report, the use of refer- 
ences 2 and 3 in their present form is justified only in cases for which 
isentropic flow may be assumed. 


According to reference 4, the velocity components at the cone 
surface are, using the notation of references 1 to 3 for quantities 
tabulated therein, 


u l(<p) 

u 


1 + a — cos <p + 
u 


UO 1 

X 

1 

/ u 2 

1 

X 

l\ 

u + 2 

— cot 0 
u 

+ 2 + 


" 2 

— cot 
u 

0 + gjcos 2<P 


( 13 ) 
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w l(<P) 


u 


Z o 

= a sr sin cp + 
u 


*2 z 

= — esc © - -=b- cot 0 
u u 


sin 2q> + 


(14) 


The pressure and density are 

= l + dl COS f + 

p p 


a 


PO 


+ \ M 2 + ^ cot © + + — 

p 2 2 p Ip 2 


M 


- -=L cot cos 2<p 


(15) 


Pl(<P) 


1 + a cos <p + 


JPO 1-2 1 £ /P 2 1 

cc I — + 2 M + 2 — cot 0 + ( — + g 


M 2 


©) 


1 € 

2 ~ cot ©J cos 2<p 


(16) 


The barred quantities are those pertaining to the case of zero angle of 
attack. From equations (lOb), (l3) , and (14), 



The plus or minus sign refers to <? = 0 or n, respectively. From con- 
ventions (.4 ) , equations (l3), (l5 ) , and (16 ) , and the result of refer- 
ence 1 that l/T = u 2 /(l - u 2 ) = — g — ) 


_ 1 _ 

I-. 


T M 2 



i 

p x p 


—2 


u 


—2 
1 - u 



1 




P 


r-i 



+ 
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or, approximately, 



In figure 5 are shown k and l/T-]_ as functions of a for a cone 
of semi vertex angle of 7.5° and a stream Mach number of 3.1. From the 
tabulations of reference 2, it may be inferred that under these condi- 
tions the isentropic assumption leads to errors of less than 1 percent 
in quantities proportional to the angle of attack. 


Solut ion of Equat ions at <p = 0 

Equations (ll) have been solved, subject to boundary conditions (12), 
at cp = 0 and various angles of attack for a cone of semivertex angle of 
7.5° and a stream Mach number of 3.1 for which the values of k an d 
l/Tq are given in figure 5. The computations were carried out by 
Dr. Lynn Albers of the Lewis laboratory and are described in appendix B. 
The resulting boundary-layer profiles of meridional velocity u and 
gradient of circumferential velocity 3w/Scp are shown in figure 6. The 
curves for a = 0 are obtained from reference 7 . The profiles show 
clearly that, as the angle of attack is increased, the boundary layer 
becomes thinner on the bottom of the cone, and the shear stress at the 
wall increases . 

Skin friction. - The meridional and circumferential components of 
the viscous shear stress at the bottom of the cone surface may be written 
in coefficient form as follows: 
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where the quantities on the right are in dimensional form. Application 
of equations (id), (2), (3), (l), (5), and (lOa) yields 



(19a) 

(19b) 


(19c) 


where 

p^x 

^1 

Variation of these skin-friction coefficients with angle of attack is 
shown in figure 7 for a particular case. 


Displacement thickness . - In reference 11 it is shown that the dis- 
placement thickness A for a cone at angle of attack is the solution of 
the equation 



At <p = 0, w-^ = 0 and equation (20 ) may he solved directly 


A = 


S x + kSqj 

1 + k 


( 22 ) 


where k is defined in equation (lOb). From equations (5) and (9), with 
5 X defined in terms of a Reynolds number. 
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(23a) 


Using equations (10) and applying 1 'Hospital's rule to evaluate at 
<p = 0 the limit of the ratio w/w]_ appearing in equation (21 ) yield 



Figure 8 shows the variation of displacement thickness with angle of 
attack for a particular case and again illustrates the expected progres- 
sive shift of the boundary layer from bottom to top as the angle of 
attack is decreased. 


Failure of Method at <j> = it 

Except for quite small angles of attack, equations (ll) cannot be 
solved at the top of the cone (<p = it ) . Over part of the range of angle 
of attack, the solutions are indeterminate j and, beyond a certain angle 
of attack, the solutions do not exist at all. These properties of 
equations (ll) will be demonstrated and discussed in the following 
paragraphs . 

Asymptotic forms of equations . - The difficulties just mentioned ma y 
best be inferred from the asymptotic forms of equations (ll) at large X. 
From equation (l2a) it is clear that, for large X, f and \|r may be 
written as follows: 


f = X + F(X) 

t = *- + ? (x) 
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where F’(“) = \|r'(«) = 0. Substitution into equations 
asymptotic forms for large X 

(l + k)XF’ ' + 2F' " =0 

(1 + k)\I ' 1 + 2Y” ’ - 2^k + = § F’|l + ^ 


(ll) yields the 



Indeterminate solutions . - Consideration will now he given to the 
problem of obtaining the complementary solution of equation (24b). 
Defining a new dependent variable 

1+k ,2 

G h e 8 r 


yields the equation 

G ” - [i k + h + (^t 1 ) x2 ] g = 0 C 25 ) 

This is essentially Weber’s equation (reference 12, paragraph 16.5), and 
the asymptotic solutions are 



Thus, the asymptotic complementary solutions of equation (24b) are 



(26a) 


(26b) 


Because it is required that Y ’ (") = 0, solution (26b) is rejected 
when k > - l/3. If k < - l/3, both solutions may be retained, and an 
additional undetermined constant appears . When k -> 0 (at the bottom of 
the cone), the complete solution of equations (ll) exists and is unique, 
and solution (26b) is to be rejected in forming the asymptotic solution. 
Therefore, if solution (26b) and the associated constant must be retained 
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when k < - l/3, it is clear that the complete solution for k < - l/3 
cannot he unique. This indeterminacy has been verified by numerical 
integration of equations (ll) , as described in appendix B. 

This indeterminacy arises because essential information has been 
lost by specializing the equations to apply only in the plane of symmetry. 
When the equations are so written, it is implied .that the boundary-layer 
development in the plane of symmetry is affected only by conditions in 
that plane. The lateral region of influence of points on the body in 
the plane of symmetry grows parabolically (the shaded regions in 
sketches (a) and (b) ) , according to the law of molecular diffusion?. 



Top view of 


cone 



(<p = it) at positive angle of attack. 


when, as in the present instance, there is no pressure gradient in the 
stream direction. Fluid enters the boundary layer from the outer stream. 
If the fluid then moves laterally (because of the angle of attack) out of 
the region of influence of the plane of symmetry, as shown in sketch (a), 
the flow is uniquely determined by outer stream conditions in the plane 
of symmetry. Clearly, this is the case when <p = jt and a is small, 
and when <p = 0 and a has any positive value. When <p = jt, except 
for small angles of attack, the lateral motion of the fluid is inward 
relative to the region of influence, as shown in sketch (b). This fluid 
then brings into the region of influence of the plane of symmetry inform- 
ation concerning boundary-layer development as it proceeds around the 
cone from the bottom. Consequently, outer stream conditions at <p = jt 
may not uniquely determine the boundary -layer characteristics . at <p = it, 
and indeterminate solutions of equations (ll) may be anticipated. 

The condition k * - l/3 specifies an angle of attack such that, at 
(p = it, the outer flow streamlines are just tangent to the parabolic region 
of influence. This statement may be proved as follows: From definition 
(10b), if k = - 1/3, then w n /u-j = - 0/2. Wear (p = it (see sketch (c)). 
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V 1 = w lq> = “ Thus, when k = -l/3, w^/uq = t/2. From this 

last equation and the geometrical properties of parabolas, it may be 
inferred that the streamlines are parabolas with focus at the cone apex. 
Therefore, the situation shown in sketch (b) applies if k < - l/3 and 
equations (11) have no unique solution. 

If k > - l/3, it might be inferred that equations (ll) may be 
solved uniquely. However, the condition k i - l/3 is not sufficient, 
but, rather, only necessary for uniqueness . Inspection of the profiles 
shown in figure 6 for a = k = 0 indicates that, for 0 > k > - l/3, 
the streamlines within the boundary layer may be expected to incline more 
sharply toward the plane cp = jt than do the streamlines near the edge 
of the boundary layer, and thus may bring information from beneath the 
cone even though the outer ones do not. Therefore, the necessary condi- 
tion for uniqueness would be (see equations (2) and (lOa)): 



or, 



(27a) 


Figure 6 indicates that perhaps the maximum value of \|r '/f ' is to be 
found at X = 0, in which case criterion (27a) would become 


k > 


1 f ' ’(0) 
3 i|r"(0) 


(27b) 


Nonexistent solutions . - The solution of equation (24a) is 

1+k, 


F ! * = e 


4 




The requirement that F ! ' (*°) =0 is met only if k > - 1. If k <.- l, 
no solution of equation (24a) exists which satisfies the boundary condi- 
tions; and; therefore, the Prandtl boundary-layer equations fail to des- 
cribe the flow. This was first pointed out by Hayes in reference 5. The 
Prandtl equations differ from the exact equations essentially in that a 
thin boundary layer is assumed. Thus, if k < - 1,' the boundary layer 
cannot be regarded as thin. It may be noted that equation (22) implies 
that, as k -* - 1, the displacement thickness approaches infinity. 
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Any boundary layer grows by the entrainment of fluid at its outer 
edge. That is, fluid particles acquire vorticity by entry into the 
boundary layer. The reverse process cannot occur - fluid particles can- 
not leave the boundary layer, thus losing their vorticity. In the case 
under consideration, it will be shown that when k < - 1, the streamlines 
at the outer edge of the boundary layer would proceed outward relative to 
the boundary layer, if the boundary layer were to retain parabolic simil- 
arity. Because such a situation is physically impossible, the Prandtl 
equations fail to yield a solution. 


From reference 11, the normal velocity at the outer edge of the 
boundary layer at <P = jt is 


v l“ U 1 ( h “ A ) 



y=o 


(28) 


where h is a somewhat arbitrary definition of the outer edge of the 


boundary layer and 


fbvA 


is obtained f rom, analys is of the outer 


nonviscous flow. The equation of continuity for the outer flow, eval- 
uated at the surface of the cone in the plane of symmetry, may be 
written 



1 ^1 

^p~ + 



(29) 


With equations (28) and (29) combined, the flow inclination at the outer 
edge of the boundary layer is 


li = ^ _ h - A 

u^_ x 



(30) 


The Prandtl (thin) boundary layer may be expected to exist only if 


dh V 1 


or, with equation (30) introduced, if 


•1 


1 Sw l 
8u-l d<p 


> - 


i + 4 0* - a) 


(31) 


If parabolic similarity is assumed (h and A each proportional to 
->/x) , inequality (31) becomes 



MCA TN 2844 


17 


1 ^1 3 

0u^ dtp > 2 


or, from equation (lOb), 


k > - 1 


SEPARATION 

The critical condition k = - 1 might be expected to be of physical 
as well as analytical significance because some sort of catastrophic 
thickening of the boundary layer is implied. In particular, this critical 
condition may reasonably be supposed to be connected with the phenomenon 
of laminar separation. In order to explore this possibility, it is first 
necessary to describe qualitatively what is meant by separation in three- 
dimensional boundary -layer flow. Difficulty has been encountered in 
establishing a satsifactory qualitative criterion for three-dimensional 
separation (see, for example, reference 5) . Therefore, in the subsequent 
paragraphs, the general problem of three-dimensional separation will be 
discussed, and then the particular case of the cone at angle of attack 
will be considered. 


General Considerations 



Plane flow . - In plane flow, separation is customarily identified 
by the appearance of reverse flow 
(sketch (d)). In order to gener- 
alize this concept to three- 
dimensional flow, it is necessary 
to consider the separated region 
as a whole . In plane flow, the 
separation point of sketch (a) 
might be regarded as the forward 
boundary of a vortex sheet 
embedded, or encapsulated, with- 
in a region bounded by the body 
and a stream surface meeting the 
body (sketch (e). Sketch (e) 
shows separation followed by 
reattachment. Of course, the 
sort of separation of greatest 
engineering importance occurs 
when such an embedded vortex 
sheet rolls up to form a large 
concentrated vortex, or is shed 
as a vortex street, with the consequence that the outer .flow is greatly 
disturbed and a large pressure effect (form drag) occurs. 


jT r TTTTTTT, 
•Separation point 
(d) Velocity profiles in plane flow. 
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If a boundary-layer solution of the type shown in sketch (e) were 
obtained, there would he little theoretical expectation of "rolling up". 
That is, a vortex sheet completely embedded in a thin boundary layer 
would be constrained to remain flat against the body. However, if the 
theoretically predicted vortex sheet extends downstream into a region 
where the thin boundary-layer equations do not apply (that is, where the 
solution "blows up" predicting an infinitely thick boundary layer), then 
in that region the vortex sheet would not be constrained to lie flat and 
the rolling -up process could occur. 

The foregoing discussion seems to provide the proper description of 
what occurs behind a bluff body: A 

complete solution- of the thin laminar 
boundary -layer equations for the flow 
over a cylinder would probably yield 
a streamline pattern of the form 
shown in sketch (f). The boundary 
layer would be predicted to gain 
mass flow by entrainment as it pro- 
ceeds around the body, until it 
reaches the vicinity of the rear stag- 
nation point. There, the mass flow 
contained in the boundary layer must 
finally leave the body and proceed 
downstream. The boundary layer 
therefore cannot remain theoretically 
thin, but rather must approach infin- 
ite thickness in violation of the Prandtl assumptions . In this region, 
then, the aft boundaries of the pair of vortex sheets are free to roll 
up into concentrated vortices, thus distorting the outer flow in such a 
way that the rolling-up process engulfs most of the region which would 
otherwise be occupied by a flat vortex sheet. The leading edge of the 
sheet, however, is still constrained to lie flat against the body. 

Three-dimensional flow . - The 
at ion may be generalized to three- 
dimensional flow as follows: A 

separated region on a three- 
dimensional body consists of a 
vortex sheet embedded between 
the body surface and a stream 
surface attached to the body 
in a closed curve, as shown in 
sketch (g), which is a view of 
the body from above . The arrows 
indicate the direction of resul- 
tant shear stress at the surface 
and outside the separated region. 

The situation shown in the sketch 
would correspond to separation 


foregoing description of plane separ- 



(g) Separation in three-dimensional 
flow with reattachment * 
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and reattachment in plane flow. If somewhere within or at the boundary 
of the separated region the theoretical boundary -layer solution would 
a priori be expected to blow up, then the vortex system within the 
separated region is free to roll up into a more or less vigorous system 
of vortices. 

Thus, thin-boundary-layer theory can be used to obtain the follow- 
ing information concerning laminar separation: 

(1) The solution may establish the existence of a vortex sheet 
which is embedded in a flat bubble on the surface and which could adhere 
to the surface and remain part of a thin boundary layer, provided that 
the Prandtl equations are valid everywhere in the separated region. 

(2) The solution might predict the boundary layer to go to infin ite 
thickness somewhere in the separated region. If this happens, then the 
separated region is free to roll up, thus providing a vigorous wake 
(which, of course, is not amenable to boundary-layer theory) . 

For flow about a plane body, the boundary-layer solution is not 
needed for predicting the breakdown of the Prandtl assumptions. Physical 
considerations suffice to establish where (at stagnation point of outer 
flow) and when (always) the breakdown occurs. In three-dimensional flow 
this is not always so clear, and at least certain features of the solu- 
tion are required to be known. In order that the boundary-layer equa- 
tions be applicable, the solution must be such that the boundary layer 
entrains fluid (that is, flow streamlines enter, but do not leave, the 
boundary layer at its outer edge). In a cartesian system, where h(x,z) 
is the outer edge of the boundary layer, this requirement may be written 
as 


Equations (32) and (28) 
in the form 


Ufhx + w 1 h z > (vq) 


y=h 


(32) 


may be combined with the equation of continuity 



to yield 



or, in vector notation. 



div 


(33) 
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where A is the displacement thickness, and qp is the velocity vector 
in the outer flow evaluated at the body surface. 


In many cases , circumstances may be found for which inequality (33) 
cannot be satisfied. For example, for plane incompressible flow about a 
cylinder and, as is customary, with u and x defined parallel to the 
surface, inequality (33) may be written 


a 

dx 


(h - A) 

h - A 


> 


duq 


dx 

U 1 


As the rear stagnation point of the outer flow is approached, up tends 
to zero, while -dup/dx remains finite and positive. Therefore, since 
h - A by definition must be greater than zero, d(h - A)/dx must 
approach infinity in clear violation of the boundary-layer assumptions. 


Separation on Cone at Angle of Attack 

In the previous discussion, it was concluded that separation involv- 
ing a strong vortex pattern occurs if a tentative boundary-layer solution 
predicts an embedded vortex sheet coupled with a local breakdown of the 
assumption of a thin boundary layer. In the case of the cone, 
inequality (33) may be used to predict the circumstances under which the 
boundary layer may not be regarded as thin: When the fact of parabolic 

similarity is introduced (h - A proportional to V*)> inequality (33) 
■becomes 


a. (h - a) 

h - A > 


-Wn 



W 1 


(34) 


This inequality indicates infinite boundary-layer thickness only when 
Wq = 0 and, even then, only if Wq is negative (which is true at the 

top of the cone, <p = it) and larger in magnitude than ^ 0uq . This is 

true only for angles of attack larger than that for which k = -1 (by 
equation (10b)). When the angle of attack is smaller than this critical 
value, the right member of inequality (34) is always negative, a-nd 
<3(h - A)/c Sep may be considered to vanish by symmetry at <p = it without 
violating inequality (34) . 

The foregoing result may be explained bn physical grounds as follows: 
As the boundary layer proceeds around the cone, it entrains fluid which it 
then conveys towards the top (symmetrically, from both sides of the cone). 
In the plane cylinder case, the fluid similarly conveyed must finally erupt 
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from the boundary layer -when the stagnation point is reached. However, 
on the cone, the boundary layer grows parabolically along generators? 
and, hence, if the cross-flow is not too strong (small angle of attack), 
the fluid brought to the top may simply become part of the growing 
boundary layer. For larger angles of attack, the boundary layer cannot 
grow at a rate sufficient to absorb the additional fluid, and eruption 
occurs with the consequent breakdown of the thin-boundary-layer 
assumptions . 

Accordingly, it is proposed that when the angle of attack is less 
than that for which k = - 1, a thin boundary layer covers the cone 
(fig. 3(a)). For larger angles of attack, any vortex sheet present will 
roll up to form attached lobes (fig. 3(b) ); for still larger angles of 
attack, a vortex street is produced. 

Thus when k < - 1 (angle of attack greater than that for which 
k = l), strong viscous cross forces (viscous lift) on the cone may be 
expected. These forces are discussed by Allen in reference 13. Of 
course, a weakness of the present analysis is that no indication is 
given as to the strength of the rolled-up vortex system because, when 
k = - 1, the presence of an embedded vortex sheet over the top part of 
the cone has not been established. It seems likely that such a vortex 
sheet does exist because a rather strong adverse pressure gradient 
(fig. 2(b)) always exists when k = - 1. In fact, it may be shown (most 
conveniently by evaluating equation (7a) at X = •) that an adverse pres- 
sure gradient exists when k < - 2/3. 

Equations (35), (13), (14), and (lOb) give the minimum angle of 
attack for which the boundary layer may be expected to break away from 
the surface and to form lobes at the top of the cone. This angle is 
presented as a function of stream Mach number and vertex angle. Fig- 
ure 9 shows the results of such a calculation. The critical angle of 
attack is given as a ratio of angle of attack to semi vertex angle for 
convenience. The results show that, in general, separation involving 
lobes occurs later (in terms of relative angle of attack a/©) for the 
smaller vertex angles, particularly at higher Mach numbers. Figure 9 
indicates the possibility of rather profound qualitative differences in 
the flows at high Mach number about cones of different Vertex angles. 

The foregoing interpretation of the critical condition k = - 1 is 
supported by the experimental result shown in figure 4. From figure 5, 
k = - 1 when a = 6.2°, under the conditions of the test. Figure 4(b) 
shows the measured total head rising as the angle of attack is increased 
beyond 6.2°. Presumably this effect is caused by the the induced field 
of the symmetrical pair of vortex lobes sweeping away the thick boundary 
layer between, thus reestablishing a thin boundary layer at the top of 
the cone. 
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It may "be of interest to note that if similarity also holds for the 
turbulent boundary layer on a cone, and the similarity lav is nearly 
linear (rather than parabolic as in the laminar case), separation would 
first appear at a higher angle of attack than in laminar flow. In fact, 
equations (lOb^ and (3l) or (33) would yield the criterion k = - 4/3. 

For the boundary layer produced by supersonic flow over any smooth 
conically symmetric body in supersonic flow (such as a cone of elliptic 
cross section), inequality (33) and the condition of parabolic boundary- 
layer similarity may be used to find a criterion equivalent to k = - 1 
for the maximum angle of attack consistent with a thin boundary layer. 


CONCLUSIONS 

The laminar boundary-layer flow about a circular cone at large angles . 
of attack to a supersonic stream has been analyzed in the plane of symmetry 
with the following results: 

1. At the bottom of the cone, profiles of meridional velocity and of 
the gradient of circumferential velocity were determined and showed the 
expected tendency of the boundary layer to become thinner on the underside 
of the cone as the angle of attack is increased. 

2. At the top of the cone, except for very small angles of attack, 
the analysis (which is restricted to the plane of symmetry) failed for 
the following reasons: 

(a) For angles of attack greater than some rather small value, the 
boundary layer brings information from beneath the cone into the vicinity 
of the plane of symmetry at the top. Therefore, the analysis, which deals 
only with the plane of symmetry, yielded indeterminate solutions. 

(b) For angles of attack greater than some angle (roughly of the 
order of the cone semivertex angle), no boundary -layer solution is pos- 
sible. The characteristics of the outer flow and the known parabolic 
similarity of the boundary layer would together imply that, beyond this 
critical angle, there would be a component of flow leaving the boundary 
layer. This is physically impossible, since a boundary layer always 
entrains fluid. Thus, beyond the critical angle of attack, no solution 
can exist for equations which presume a thin boundary layer. 

For three-dimensional flow it is proposed that a separated region be 
regarded as a vortex sheet embedded in the boundary layer, remaining flat 
against the body if the assumption of a thin boundary layer is valid 
throughout the region. If, however, the boundary -layer assumptions break 
down anywhere in the separated region, it is inferred that the vortex 
sheet may roll up to form strong vortices which may either remain attached 
or be shed as a vortex street* 
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On the cone, therefore, the critical angle of attack beyond, which 
no boundary -layer solution is possible at the top of the cone represents 
the maximum angle of attack for which the boundary layer is everywhere 
thin or, alternatively, the minimum angle of attack for which major dis- 
ruption of the flow may be expected because of the formation of strong 
vortex lobes . Beyond this angle of attack, strong viscous cross forces 
may be anticipated. 

A similar criterion could easily be obtained for the boundary layer 
on any smooth conically symmetric body in supersonic flow. 

The assumption of a suitable similarity law suffices to establish 
a similar criterion if the boundary layer on a conical body is turbulent. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 15, 1952 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 

C constant appearing in temperature -viscosity relation 

(equation (5a)) 

Cf component of skin-friction coefficient in x-direction 

X x 

component of skin-friction coefficient in ^-direction 
<P 

Cp specific heat at constant pressure 

F(X) function appearing in asymptotic representation of f 

(equations (2b) 

f(\><p) function related to meridional velocity u by equation (2) 

function related to circumferential velocity w by 
equation (2) 

h height of outer edge of boundary layer 

k related to circumferential gradient of circumferential 

velocity in plane of symmetry (equation (lOb)) 

M Mach number 

p static pressure 

velocity vector at outer edge of boundary layer 

R x Reynolds number, PqUqx/ \ij_ 

R/^Rg^Rs Reynolds numbers, Pl u l A /Mq^ Pl u l 5 x/^1> Pl u l^p/^ 

^ respectively 

T absolute static temperature 

u meridional component of velocity 

v component of velocity normal to surface 

v circumferential velocity component 
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x coordinate along generators of cone 

y coordinate normal to surface 

a angle of attack (positive as shown in fig. l) 

T ratio of specific heats 

A displacement surface 

& x mass-flow defect associated with meridional velocity pro- 

file (equation (2l)) 

&(p mass -flow defect associated with circumferential profile 

(equation (2l)) 

0 semivertex angle of cone 

9 sine of semi vertex angle of cone 

X dimensionless variable (equation (3)) 

p coefficient of viscosity 

p density 

<p angular coordinate around cone 

f(X) function appearing in asymptotic representation of i|r 

(equat ions (24 ) ) 

t(x) function related to circumferential velocity w in plane 

of symmetry by equation (lOa) 


Subscripts: 

max maximum 

r reference condition, nonviscous flow at surface, at <p = 0 

or cp = jt, whichever is appropriate 

1 evaluation at outer edge of boundary layer (alternatively, 

nonviscous flow at surface) 

Subscript notation for partial differentiation has been used 
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Superscripts: 

' Primes denote ordinary differentiation with respect to \ or cp 

— Bar over quantity indicates evaluation of nonviscous 

flow at surface when cone is at zero angle of attack 
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NUMERICAL SOLUTION OF DIFFERENTIAL EQUATIONS 
By Lynn Albers 

The two simultaneous nonlinear ordinary differential equations (lla) 
and (lib) together with boundary conditions (12) constitute a two-point 
boundary value problem. The method of numerical solution used applies 
directly only to problems fow which all boundary conditions are specified 
at a single initial point (the origin, in the present case). Each numer- 
ical integration was therefore performed starting with boundary condi- 
tions (l2a) and (12b) and a tentative specification of f''(0) and 
t''(o). In each case, such integration was carried out for a sufficient 
variety of conditions f''(0) and ' ' (o) so that the correct set of 
initial conditions yielding the proper behavior at X = ~ (boundary con- 
dition (12a)) could be inferred to the desired degree of accuracy. 

Integration was performed according to the following basic scheme: 
With the values of f'''(X) and T|r'''(x) given at five closely spaced 
values of X, fourth-degree polynomials may be passed through the two 
sets of values of f' 11 and Then, if f , f ' , f 1 ' , ijr, \|r ' , and 

t' ' are known at the fifth point, the polynomial representations of f ' ' ' 
and maybe integrated to yield f, f', f' 1 , ijr, \|r', and * ' at 

the next (sixth) point. These quantities may then be substituted into 
differential equations (ll) to yield f * 1 1 and V ! ' at the sixth point. 
In this way, the solution may be extended one step at a time, in each 
step by use of the solution at the five previous points. In order to 
begin this procedure, the solution must first be found at five points 
starting at the origin and must be subject to boundary conditions (l2b) 
and (l2c) and the tentative selection of f'(o) and ^''(o). 

This preliminary calculation was done in the following manner: 
f' ' ' (0 ) was calculated directly from equations (ll) and (12) and was 
used as an initial estimate of f 1 ' ' at the next four points. Given 
f(o), f'(o), and f''(o), the values of f, f', and f' 1 were computed 
at the second point by integrating a fourth-degree polynomial passed 
through the five values of f ' ' '. In a similar manner, ijr, i)f ' , and \|/' ' 
at the second point were found. Direct substitution into equations (ll) 
then yields improved estimates for f ' ' ' and \|/' ' ' at the second point 
and thus an improved polynomial representation of these functions which 
may be used to obtain values of f, f', f ' \jr, ilr', and \jr' ' at the 
second point, and so forth, until improved values have been obtained at 
the fifth point. This procedure was repeated in an iterative manner 
until convergence was obtained at each of the five initial points. 

All calculations were performed on the IBM Card Programmed Electronic 
Calculator. Results are considered correct to four significant figures. 
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